Magnetic Properties of the Intermediate State in Small Type-I Superconductors 
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We present simulations of the intermediate state of type-I superconducting films solving the time 
dependent Ginzburg-Landau equations, which include the demagnetizing fields via the Biot-Savart 
law. For small square samples we find that, when slowly increasing the applied magnetic field 
Ha, there is a saw-tooth behavior of the magnetization and very geometric patterns, due to the 
infiuence of surface barriers; while when slowly decreasing Ha, there is a positive magnetization and 
symmetry-breaking structures. When random initial conditions are considered, we obtain droplet 
and laberynthine striped patterns, depending on Ha- 

PACS numbers: 74.25.-q,74.25.Ha,75.60.-d 



In 1937 Landau modeled the intermediate state (IS) 
in thin slabs of type-I superconductors, assuming a pe- 
riodic structure of alternating stripes of normal and su- 
perconducting phases.^ Direct experimental observation 
of the IS revealed that, while in some cases its structure 
resembled the Landau picture, very complex patterns 
and history dependence were usually seeni^*^*^ Similar 
type of complex structures were later observed in two- 
dimensional (2D) systems where there is a competition 
among interfacial tension and long-range interactions^ 
like thin magnetic films, ferromagnetic fluids, Langmuir 
and lipid monoloyers, and self-assembled atoms on solid 
surfaces.^ Labyrinthine patterns, and a transition from 
structures of droplets to stripes are typically observedi^ 
The rich physics found in these 2D systems has moti- 
vated a renewed interest in the study of the IS in type-I 
superconductors in several recent experimentsi2i2ii2iii 

Most of the theoretical progresa^^i"d^i^^d^i^'^d^i^^i^° 
has been made by modeling the IS with periodic ar- 
rays of normal and superconducting phases. Recently, 
a current-loop modeU^ which allows to describe simple 
non periodic patterns has been introduced, but a fully 
consistent theoretical descriptioniS of the experimental 
patterns is still needed. Another important problem 
not addressed neither experimentally nor theoretically in 
type-I superconductors is the study of the IS in samples 
of sizes comparable with the expected periodicity of the 
patterns, while interesting "mesoscopic" behaviors have 
been found in type-II superconductors with sizes of the 
order of few times the magnetic size (A) of vorticesiSi 

In this paper we report detailed simulations of small 
square type-I superconductors by solving the time de- 
pendent Ginzburg-Landau (TDGL) equations for slabs 
of thickness d. We consider the "non-branching case" , 
where d <^ ds ^ 800(^ — A)^^ which can be well approx- 
imated by reducing th e equa tions to a 2D problem, as 
done for example in Ref fl^flSl We therefore assume that 
the current density J and the order parameter 5" can be 
replaced by their average over z for —d/2 < z < d/2, i.e. 
J(R,z) J(R) and *(R, z) -> *(R), with R = {x,y) 
the in-plane coordinate. The TDGL equations ,J'^^ in 



the gauge where the electrostatic potential is zero, are 
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Here V, A, J are 2D in-plane vectors, tq = ^ lO, D 
is the diffusion constant, cr„ the normal state conductiv- 
ity, A(r) the penetration depth and ^(T) the coherence 
length. These 2D approximated TDGL equations cou- 
ple with the perpendicular component of B. Following 
Refi22- we express the z-averaged sheet current J by a 
scalar function g: J(R) = V x zgCR). This guarantees 
that V- J = 0, the physical meaning of g(R) being the lo- 
cal magnetization or density of tiny current loops. Next 
one relates ^(R) with Bz{r) = V x A\z at z = and the 
applied field Ha by means of the Biot-Savart lawi2^ 

Bz(R,z^O)^Ha + ^J g(R, R')g(R')'^'R' (3) 

The kernel Q satisfies (5(Ri,R2) = Q(R2 — Ri); 
Q{R > d) = -d/R^; and / d'^IiQ{K) = (due to 
flux conservation). To a good approximation the ker- 
nel can be given by (see K. Maki in Ref. 12): Q(R) ~ 
47r(5(R) - d/[|R|2 -I- d'^/4f^'^. The boundary conditions 
are (V - ij^A)\j_^ = and g\b = 0. We solve the 
equations with a finite difference scheme with discretiza- 
tion Ax = Ay = 0.5^(0), using link variables to main- 
tain the gauge invariancei^i We normalize time by ta = 
AnanX'^/c^, A by 7?c2(0)C(0) , T (temperature) by T^, 
and take Tc/tQ — 12. To obtain g{R) we invert (3) using 
the conjugate gradient method as done in R,efl23l 

We show results for square samples of size L x L with 
L = 256C(0), thickness d = 40^(0); k ^ \/C ^ 0.6; 
T — O.STc and time step At = 0.25 (we also obtained sim- 
ilar results for L = 120^(0) -256^(0), d = 10^(0)^40^(0) 
and K — 0.4 — 0.6). We performed simulations of the 
intermediate state following three different procedures: 

(i) slowly increasing the magnetic field from Ha = 0, 

(ii) slowly decreasing the magnetic field from the normal 
state {Ha > He) and (iii) starting from random initial 
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FIG. 1: (a) -{B - Ha)/Ha{T) and (b) the number of flux 
quanta obtained increasing (closed circles), decreasing (open 
circles) the external magnetic field h = Ha/ Hc{T), and with 
random initial conditions (crosses). 



conditions for each value of Ha- The global results for 
the three cases are summarized in Fig. 1. We show the 
apparent magnetization, AuMa = {Bz) — Ha (real mag- 
netization is AnM — B — H, but Ma is what can be 
determined experimentally), in Fig. 1(a), and the number 
of flux quanta inside the sample, N^o = §{A + j^)dl, 
in Fig. 1(b), as a function of h = Ha/Hc{T). 

(i) Slowly increasing the magnetic field. We start from 
Ha ^ with a state with |*(R)P = 1 and B{R) = 
and increase Ha in small steps, after reaching a station- 
ary state for each field (when \ AE/E\ < 10~^; where AE 
is the change in energy between consecutive time steps). 



We observe a Meissner state for Ha < Hp = 0.56Hc. Sur- 
face barriers preclude the penetration of flux below Hp. 
The surface barrier in macroscopic type-I superconduc- 
tors can lead to relatively large first penetration fields Hp, 
which depend on the sample shape and dimensions. ^'^^ 
Here the smallness of our system strongly enhances this 
effect. We observe that at H > Hp four long chunks of 
the normal phase, carrying hundreds of flux quanta, en- 
ter from each side of the square and equilibrate in the 
pattern shown in Fig. 2(a). At a higher field -Hp, 2, some 
other four chunks of flux enter and form the pattern seen 
in Fig. 2(b). Further increasing the field, more complex 
structures form, as shown in Fig. 2(c) and (d). The nor- 
mal domains tend to stay in the centre of the sample, 
leaving a flux-free zone near the edge. Above H^ flux 
has fully entered inside the system and there is only 
surface superconductivity until i?c3 > He. The inter- 
nal structure of the domains is detailed in Figs. 2(e) and 
2(f) which show transversal cuts of |\['(R)|^ and Bz(R) 
taken at the center of one of the faces for h = 0.58. In 
Fig. 2(e) we see that in the normal regions there is a 
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FIG. 2: (a)-(d) Spatial patterns of |$(r)p obtained contin- 
uously increasing h — Ha/ He from h = 0. (a) h=0.58, (b) 
h=0.65, (c) h=0.72 and (d) h=0.78. Gray scale ranging from 
black for [^(r)^ = to white for |*(r)p = 1. (e) and (f) 
show transversal cuts of |^(r)p and (Bz) taken at the cen- 
ter of one of the faces for h = 0.58 (Fig. 2(a)). Parameters: 
K = 0.6, d = 40^(0) and sample size 256^(0) x 256^(0). 



sharp drop to zero of l^tp and that B^ ^ Hc{T) (see 
Fig.2(f)). The global structure of the patterns of the IS 
observed in Figs.2(a)-(d) follow the geometry of the sam- 
ple and have the symmetry of the square. In general, we 
find that the entrance of the normal phase occurs only 
for discrete values of penetration fields Hp ^ where sev- 
eral fiux quanta are nucleated at the four sides of the 
square, while for Hp i < Ha < there is no flux 

entrance. This shows up as a saw-tooth behavior in the 
magnetization in Fig. 1(a) and as a series of plateaus and 
jumps in the number of flux quanta vs. Ha in Fig. 1(b). 
This type of behavior is similar to the results observed 
in mesoscopic type-II systems^^ which also show a saw- 
tooth behavior of the magnetization. However, while in 
Ref. [21] each jump in Ma corresponds to the entrance of 
one quantum of fiux (one vortex), here at each jump in 
Ma several hundreds of flux quanta have entered. Our re- 
sults suggest the existence of a "mesoscopic-like" behav- 
ior in small type-I samples. This novel behavior appears 
when the linear size L of the sample only allows for a 
small number of normal domains inside the system. This 
means that L is not more than one order of magnitude 
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larger than the periodicity of the patterns at intermediate 
fields. Indeed, we have found similar "mesoscopic-like" 
behavior for sizes in the range L < 256^(0). 
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FIG. 3: Spatial patterns of obtained decreasing h from 

the normal state at /i 3> 1. Same parameters as in Fig. 2. 



(ii) Slowly decreasing the magnetic field. We start from 
Ha > He with a state with = and B = Ha and de- 
crease Ha in small steps, after reaching a stationary state 
for each field. The resulting intermediate state patterns 
are shown in Fig. 3. When Ha ~ He the superconduct- 
ing phase enters into the sample and the total number 
of flux quanta is reduced. At first, the superconducting 
phase forms four chunks embedded within the normal 
phase, which nearly follow the square symmetry of the 
system, as shown in Fig. 3(a). However, we observe that 
when decreasing the field the square symmetry is always 
broken in the patterns. The breaking of symmetry is 
stronger the further we decrease the field. In this way, 
labyrinthine patterns are formed at mid-range fields, as 
can be seen in Figs. 3(b) and (c). In general, the expulsion 
of fiux occurs gradually when decreasing Ha as shown in 
Fig. 1(b). For low fields we see that thin stripes of nor- 
mal phase are trapped within the sample, as shown in 
Fig. 3(d). The difficulty for expeUing fiux is due to the 
surface barrier and results in a positive magnetization as 
a function of h as shown in Fig. 1(a). Even aX h — 0, a 
small amount of flux remain trapped in the sample, as 
evidenced in Fig. 1(b), where the number of fiux quanta is 
finite at ft, = 0, and in Fig. 1(a), where Ma > at h — 0. 
In experiments in macroscopic samples it has been ob- 
served that some trapped fiux remains in the system at 
h = when decreasing the field"^ and in some Sn films a 
positive magnetization has been obtained when decreas- 
ing Ha^ It is interesting to mention that a similar pos- 
itive magnetization was observed in mesoscopic type-II 
superconductors^^ when decreasing Ha, and attributed 
to the importance of surface barriers. 



(a) h=0.33 



(c) h=0.50 




(b) h=0.42 



(d) h=0.75 



FIG. 4: |^'(r)p patterns obtained using random initial con- 
ditions. Same parameters as in Fig. 2. 



(Hi) Random initial conditions. The structures of the 
IS discussed above, obtained either increasing or decreas- 
ing the magnetic field, are strongly influenced by the sur- 
face barriers and/or the geometry of the small sample 
simulated. In a fllm of large linear size L the demagne- 
tization factor N is such that 1 — iV oc d/L w 0, and 
therefore we expect that B « Ha. To obtain stationary 
states more typical of the bulk behavior of large samples, 
we start with a initial condition with random values of 
A and such that we satisfy {B) — Ha from the start, 
and that the initial state is superconducting in average, 
(l^'P) > 0. We performed simulations with this initial 
condition for different values of h, and let the system to 
evolve for each case, using a stronger criterion for assum- 
ing stationarity: \AE/E\ < 10~^. We obtain that in 
the stationary state most of the flux remains inside the 
sample and (Bz) — Ha w as can be seen in Fig. 1(a). 
The structures obtained are shown in Fig. 4. For low 
fields, we observe in Fig. 4(a) that the intermediate state 
consists of almost circular droplets of the normal phase. 
For higher fields, the droplets start to coalesce into long 
lamellar-like domains, as seen in Fig. 4(b). At intermedi- 
ate fields, as shown in Figs. 4(c), labyrinthine patterns of 
stripes of the normal phase are formed. For high fields 
close to He we observe almost circular droplets of the 
superconducting phase embedded within the mostly nor- 
mal phase, see Fig. 4(d). Similar type of structures, with 
droplets of one or the other phase at low and high fields, 
and with labyrinthine patterns of stripes at mid range 
fields has been observed experimentally for example in 
Figs.2.8(a)-(f) of Ref0 for a lead thin film. An impor- 
tant feature we find in our small system is that there is 
a thin layer of superconducting phase at the surface (see 
Fig. 4), which allows for the screening Meissner currents. 
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FIG. 5: (a) Spectral intensities obtained from [^(r)] at 
different magnetic fields, a maximum at ko is observed. The 
inset shows the Fourier transform of Fig. 4(a). (b) Periodicity 
a — 27r/fco of each structure as a function of h. The lines 
correspond to the Landau and the Goren-Tinkham models. 



We analyze the structures obtained in Fig. 4, by cal- 
culating the spectral transform of the superconducting 
order parameter, /(k) = |/ dr|5'(r)p exp(ik ■ r)| , which 
is shown in Fig. 5 (a). The non-periodicity and complex 



structure of the patterns results in very broad maxima 
in /(k) at finite wave vectors kg — 27r/a which define 
a typical length scale a. In the case of low and large 
fields, a would correspond to the typical distance between 
droplets, while for mid-range fields, a would correspond 
to the average widths of the stripes in the labyrinthine 
patterns. We plot a in Fig. 5(b) and compare it with the 
Landau model of stripes^ and with a model of Goren 
and Tinkham for a periodic array of droplets or "flux 
spots" We see that the Landau model agrees quali- 
tatively with the results obtained at mid-range h (for 
these fields the patterns of Fig. 4 can have a mixture of 
"stripes" with a few droplets, which make a smaller than 
the Landau value). On the other hand, at low h the 
Goren-Tinkham model does not agree well with the size 
of the droplets obtained. Also in some experiment3iii£ 
it has been found a departure from the Goren-Tinkham 
model at low fields. 

In conclusion, our simulations predict that the strong 
influence of the surface barriers in small type-I samples 
will lead to a saw-tooth behavior of the magnetization 
and very geometric patterns when slowly increasing Ha, 
and to a positive magnetization and symmetry-breaking 
structures when slowly decreasing Ha- These results sug- 
gest the existence of a "mesoscopic-like" behavior in the 
IS when the sample linear size is of the order of a few 
times the periodicity of the patterns. It will be interest- 
ing if experiments on small samples of type-I supercon- 
ductors could be performed. 
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